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ABSTRACT 

In order to test if there is energy transfer between dark energy and dark mat- 
ter, we investigate cosmological constraints on two forms of nontrivial interaction 
between the dark matter sector and the sector responsible for the acceleration 
of the universe, in light of the newly revised H{z) data, the baryonic acous- 
tic oscillation (BAO) peak detected by large-scale correlation function of lumi- 
nous red galaxies from Sloan Digital Sky Survey (SDSS), the cosmic microwave 
background (CMB) observations and the most recent 557 Union2 SNe la data. 
For the 7mIDE model in which the interaction term expresses as Q = S'jrnHpm 
with the corresponding constant 7^ quantifying the extent of interaction, we 
obtain 7™ = -0.0131°:°}?, 7^ = -0.0091°;°}^ and 7^ = -0.0111°;°}? with 
H{z)+BAO+CMB, SNe la+BAO+CMB and H{z)+SNe la+BAO+CMB, re- 
spectively. For the jiJ.DE model in which the interaction term takes the form 
Q = S'jdHpx with the corresponding constant •jd quantifying the extent of 
interaction, the best-fit values are 7^ = -0.0401o;q42, 7d = -0.023lo;o4o and 
7d = — 0.030lo;g37. These results indicate that, on the one hand, the concordance 
ACDM model without any interaction remains a good fit to the recent observa- 
tional data; on the other hand, the interaction that the energy transferring from 
dark matter to dark energy is slightly favoured over the interaction delivering en- 
ergy from dark energy to dark matter, which is in agreement with the results by 
using other observational data such as the 182 Gold SNe la and 397 Constitution 
SNe la samples in the previous works. 

Subject headings: cosmological parameters - (Cosmology:) dark energy 



1. Introduction 



Many astrophys i cal and cosmological observations, such as Type la Supe rnovae (SNe 
la) jRiess et al.lll998l : IPerlmutter et al.lll999HRiess et al.lbood : iKnop et al. l l2007[ ). co smic mi- 
crowave background (CMB) from Wilkinson Microwave Anisotropy Probe (WMAP) (jSpergel et al 
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20031 . l2007t [Komatsu et aPbood 1201(1 ) , and the large scale st ructure from Sloan Digital Sky 



Survey (SDSS) (iTegmark et al.l l2004l : lEisenstein et al.l |2005| ) have indicated that the uni- 



verse is undergoing an accelerating expansion, which consequently leads to the assumption 
of the existence of dark energy, an exotic energy with negative pressure and causes an ac- 
celerating expansion of our universe at late times. The most simple candidate for these 
uniformly distributed dark energy is considered to be in the form of a cosmological con- 
stant (A) with a equation of state w = Pa/ Pa = —1. However, the corresponding ACDM 
model is always entangled with the coincidence problem: The matter density pm decreases 
with the expansion of our universe with and the density of cosmological constant pA 
does not change with the expansion of the universe, whereas the dark energy density is 



include the scalar field models with dynam ical EoS (e. 
Ratra fc Peeble^llQSsi B 



(lArmendariz-Picon et al. 



2005 



2006 



dwe. 



2001 



1 et al 



f 
1^ 



(Wetterich 


1988; 


20031). k-essence 


3del ( 


Feng et al. 



Guo et al.ll2005al )). have been proposed to alleviate the coincidence problem, the 



nature of dark energy is still unknown. 

It is natural to consider the possibility of exchanging energy between dark energy and 
dust matter through interaction term. The interacting dark energy model was first intro- 
duced to test the coincidence problem, in whic h pm could decrease with the expansion of 
our universe slower than a~^. lAmendolal ( 120001 ) investigated a coupled quintessence (CQ) 



model by assuming an exponential potential and a li near coupling, be si des many other 



backg r ound and perturbation con straints carried out by lBean et al.l (|2008[ ): iLa Vacca. et al. 



(120091 ) : lDe Bernardis et al.l ( 120111 ). On the other hand, it is always assumed that dark energy 
and dust matter exchange energy through an interaction term Q 



px + {px + px) = -Q, 

Pm + S-ffpm = Q, 



which preserves the total energy conservation equation pt. 
teraction term Q is ex 



at 



tensiyely considered in the literature ( Copeland et al. 



2001; 


Chimento et al. 


2003; 


Setare 


2004; 


Cai & Wane 


2005 


2005; 


Wei et al. 


2005 




Wei & Cai 


2006; 


Szvdlowski 


20061; 


2007; 


Chen et al 


2009. 


201C 


i). If Q is a non-zero function 



+ {Ptnt + V tnf) 



(1) 

0. The in- 



Guo et al 



2005b 



Bertolami et al. 



9981: 



Dalai et al. 



Noiiri fc Odintsov 



2007 



Guo et al. 



makes pm and pd to deviate from the standard scaling. The simplest assumption is 

Q = ^iHp, (2) 
where p is the density of du st matter prp ( Wei fc Zhang 2007al lbl: Izhang fc~Ti 2010) or the 



density of dark energy px (IZhang fc Lil |2010| ) with the corresponding constant 7.^ or 7^ 
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quantifying the extent of interaction between dust matter and dark energy. We assume 
that the EoS of dark energy wx = Px/px is a constant in spatially flat FRW metric. When 
working out the value of 7, we can see the extent of interaction and transfer direction between 
dark energy and dark matter. The ACDM model without interaction between dark energy 
and dark matter {Q = 0) is characterized by 7 = 0, while 7 7^ and denote non-standard 
cosmology. For 7 < 0, the energy is transferred from dark matter to dark energy; while 
for 7 > 0, the energy is transferred from dark energy to dark matter, and the coincidence 
problem can be alleviated. 

For the observational data, the observational Hubble parameter data H{z) have become 
an effective probe both in cosmology and astrophysics compared with the SNe la data, the 
CMB data and the baryonic acoustic oscillation (BAO) data and it is more rewarding to 
investigate the observational H{z) data directly. The reason is quite simple, it is obvious that 
these probes all use the distance scale (e.g., the luminosity distance rf^, the shift parameter 
R, or the distance parameter A) measurement to determine cosmological parameters, which 
needs the integrate of the Hubble par ameter and ther efore lose the fine structure and some 
more important information of H{z) (ILin et al.ll2009l ). The Hubble parameter depends on 
the differential age as a function of redshift z in the form 



which provides a direct measurement for H{z) through a determination of dz/dt. Ijimenez et al, 
(120031 ) demonstrated the feasibility of the method by applying it to a 2; ~ sample. 
By using the differential ages of passively evolving galaxies determined from the Gem- 



ini Deep Deep Survey (G D DSl ([Abraham et a. 



20041 ) and archival data (ITreu et al.l 12001 



2OO2I : iNolan et al.1l2003al E]). Isimon et al.l ( |2005l ) determined 9 H{z) data in the range < 
z < 1.8. These obs ervational H(z) data we r e also used to con strain the parameters of 



cosmological models (jSamus. 



works include IWei fc Zhang 



lia &: Ratra 
j2007Jbll: 



2006 



Wu fc Yu 



Kurek fc Szvdlowskil J2008fi: ISen fc Scherreil J2008f ): IXu et all (l2008h : Izhang fc Zhul 120081 ) 



(l2007a 



Yi &: Zhang 2007 ) and some other r elevan t 



azkoz fc Maierottol (120071) : 



for examples. Recently, IStern et al.l (120101 ) obtained the H{z) data at 11 different redshifts 
obtained from the differential ages of red-envelope galaxies; and other two Hubble parameter 
data at 2; = 0.24 and z = 0.43 were determined by lGaztaiiaga et al.l (120091 ) from observations 
of BAO peaks. Some recent works using these newly H(z) dat a for cosmological constrain t 



can be found inlGong et al.l (|2010r ) : iLiang:. Wu fc Zhand (|2010l): ICao. Zhu &: Liang 



Ma fc Zhanj (l201lh : IXu fc WangI (l201oh : 



Zhai et al 



(I2OI0I ): Izhang. Ma fc LanI (I2OI0I ). 



2Qlll); 



In the previous works, IWei &: Zhangl (j2007a[ ) compared the 9 observational H{z) data 
with some cosmological models with/without interaction between dark energy and dust mat- 
ter and found that the H{z) data points with fairly large errors cannot severely constrain 
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model parameters alone. In this paper, we focus on the newly H{z) data to study the inter- 
action between the dust matter and dark energy and test the cosmic coincidence problem. 
In order to break the degeneracy of model parameters, we also add the baryonic acoustic 
oscillation (BAO) peak detected by larg e-scale correlation fun ction of luminous red galaxies 
from Sloan Digital Sky Survey (SDSS) (lEisenstein et al.l 120051). the cosmic m icrowave back- 
ground (CMB) detected by the 7 - year WMAP data (IKomatsu et al.ll2010l ) and the newly 
revised Union2 SNe la data set (lAmanullah et al.ll2010l ). This paper is organized as fol- 
lows: In section [21 we introduce the observational data including the H{z), BAO, CMB and 
SNe la data. In section [3l we derive two Hubble parameters and perform a Markov Chain 
Monte Carlo analysis spanning the full parameter space of the model using different data 
sets. Finally, we summarize the main conclusions in Section HI 



Observational data 



In this section we will list the cosmological observations used in our calculations: H{z), 
BAO, CMB as well as t he SNe la observa tions. We adopt the H{z) data at 11 different 
redshifts obtained in Ref. Istern et al.l J201oh . an d two H(z) data (H(z = 0.24) = 76.69±2.32, 
and H{z = 0.43) = 86.45 ±3.27) determined by lGaztanaga et al.l (120091 ). The corresponding 
can be defined as 



Xh 



13 



a 



(4) 



hi 



where ahi is the la uncertainty in the H{z) data. 



As is known, the baryonic oscillations at recombination are expected to leave baryonic 
acoustic oscillations (BAO) in the power spectrum of galaxies. The expected BAO scale 
depends on the scale of the sound horizon at recombination, and on transverse and radial 
scales at the mean redshift zbao = 0.35 of galaxies in the survey. lEisenstein et al.l (120051 ) 
measured the quantity 



A 



EizBAoY^' 



r 1 


r^^o dzl 1 


-Zbao J 


E{z')\ 



2/3 



(5) 



The SDSS BAO measurement (lEisenstein et al.ll2005f l gives Aobs = 0.469(n,/0.98) - 0.35 ± 
0.017 where the scala r spectral index is taken to be Ug = 0.963 as measured by WMAP7 
( iKomatsu et al.ll2010[ ). In this case, x^can be defined as 



Xbao 



{A -A 



obs j 



cr. 



(6) 
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Meanwhile, the locations of peaks in the CMB temperature power spectrum in 1-space 
depend on the comoving scale of the sound horizon at recombination, and the angular 
distance to recon i bination. This is summa. rized by the so-called CMB shift parameter R 
( iBond et al.l 119971 : IWang &: Mukherjed 120061 ) which is related to the cosmology by 



R = ^/n 



mO 



dz 

W) 



(7) 



where z^ec ~ 1091.3 ( iKomatsu et al.l l2010l ) is the reds hift of recombination . The 7-year 
WMAP data gives a shift parameter R = 1.725 ± 0.018 ( IKomatsu et al.ll2010l ). In this case, 

X 



^ can be defined as 



XCMB 



{R-R, 



obs J 



(8) 



Notice that both A and R are independent of Hq, thus these quantities can provide robust 
constraint as complement to H{z) on dark energy models. 

It is commonly believed that SNe la all have the same intrinsic luminosity, and thus can 
be used as "standard candles." Recently, the Supernova Cosmology Proje ct (SCP) collabo- 
ratio n have released their Union2 compilation which consists of 557 SNe la (lAmanuUah et al. 
2010l ). The Union2 compilation is the largest published and spectroscopically confirmed SNe 
la sample to date, which are used in this paper. In the calculation of the likelihood from 



SNe la we have marginalized over the nuisance parameter (jPi Pietro fc Claeskensll2003l ): 



XsN 



C 



(9) 



where A = E ' 'M' ,B = Et l^''^'^ - fi'"' / ,C = Ei l/^l f^"""'" is the 

distance modulus obtained from observations and ai is the total uncertainty of SNe la data. 



3. Constraint on the IDE models 

The model parameters are determined by applying the maximum likelihood method 
of fit by using the Markov Chain Monte Carlo (MCMC) method. We minimize to 



deter mine the best-fit parameters and our method is based on cosmoMC (iLewis &: Bridle 



20021 ). Basically, The model parameters are determined by minimizing 



— Xh + Xbao + XCMB + XsN- 



(10) 
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3.1. The 7^ IDE Model 

If the interaction term is Q = S'jmHpm, in spatially flat FRW metric, for the the 7^ 
IDE model with a constant EoS of dark energy wx, the Friedmann equation is 



E\z) = " (1 + zY^'-^-) + 1 - " (1 + zY^'^""^^ . (11) 



The joint confldence regions in wx-lm plane with different observational data sets {H{z), 
H{z)+BAO+CMB, SNe la+BAO+CMB, and H{z)+SNe la+BAO+CMB) for the 7^ IDE 
model model are showed in Fig. 1. We also present the best-fit values of parameters with 1-a 
and 2-cT uncertainties in Tabledl With the H{z) data only (Fig. la), the best-fit values of the 
parameters (wx,7m) are wx = -2.79 and 7^ = 0.22. With if(2)+BA0+CMB (Fig. lb), 
the best-fit values at 1-a are wx = —I-IO^oIt, 7m = — 0.0131;q;o}'^. For comparison, fitting 
results from the joint data with SNe la-|-BAO+CMB are given in Fig. Ic. with the best-fit 
values Wx = and 7^ = — 0.009]';o;q}2- In Fig- Id, we show the fitting results from 

the joint data with H{z)+SNe la+BAO+CMB, with the best-fit values wx = -l.05to\l 
and 7m = -0.011^:°:°}?. 

It is obvious that H{z) only gives a relatively weak constraint on all of the relevant model 
parameters. We find that the H{z) data, when combined to CMB and BAO observations, can 
give more stringent constraints on this phenomenological interacting scenario . Moreover, the 
special case (7^ = 0, wx = — 1, corresponding to the ACDM with no interaction) is excluded 
at la for if(2r)+BA0+CMB. However, obviously, with the two joint analyses including the 
SNe la data, ACDM is still included within la error region. Moreover, comparing Fig. lb to 
Fig. Ic, we can find the confidence regions of if(2;)+BA0+CMB data are in good agreement 
with that of SNe la+BAO+CMB data; this situat ion has also been noted when constraining 



without dark energy (ICao. Zhu fc Liang II2OIII ). 



on the ACDM, XCD M scenario (IZhai et al.l 120101) and the interacting dark matter models 



3.2. The 7d IDE Model 

If the interaction term is Q = 3'~fdHpx, in spatially fiat FRW metric, for the the 7^ IDE 
model with a constant EoS of dark energy wx, the Friedmann equation is 

E\z) = (i-a^o)(i + ^)'^'+"''+'""^ 

^;^l]^0 + 7, + 7,(r]^0-l)(l + ^)3(7.+^x) 

Wx+ld ^ ' 
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Fig. 1. — The 68.3% and 95.4% confidence level contours for wx versus 7^ with the H{z) 
data (a); H{z)+BkO+CMB (b); SNela+BAO+CMB (c); if(2)+SNeIa+BA0+CMB {d). 
The dashed lines represent 7m = and wx = — 1- 
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The joint confidence regions in wx-Jd plane with different observational data sets {H{z), 
H{z)+BAO+CMB, SNe la+BAO+CMB, and if(z)+SNe la+BAO+CMB) for the 7^ IDE 
model are showed in Fig. 2. We also present the best-fit values of parameters with 1-a and 
2-0" uncertainties for the 7^ IDE model in Table [TJ With the H{z) data only (Fig. 2a), 
the best-fit values of the parameters (wxyja) are wx = —1.30 and 'jd = —0.057. With 
i7(z)+BA0-|-CMB (Fig. 2b), the best-fit values at 1-a confidence level are (wx,7d) = 
(-1.24+°;}^, -OMOtoMl)- Fitting results from the joint data with SNe la+BAO+CMB are 
given in Fig. 2c with the best-fit values (wx,7d) = (-1.02+°;°^, -0.023+°;°^^). Though the 
special case of the ACDM with no interaction is excluded at Icr for if(z)+BAO+CMB. 
Fitting results from the joint data with H{z)+SNe la+BAO+CMB given in Fig. 2d with 
the best-fit values (wxyld) = (—1-10 
within la error region. 



hO.13 
-0.13> 



0.030lj^:[!i^^) indicate that ACDM is still included 



From Fig. 1-2 and Table[T], it is obvious that the constraints on both interacting scenarios 
favor 7 < 0, which indicates that the energy is transferred from dark matter to dark energy 
and the coincidence problem is quite severe, a result consistent with the previous results by 



(Chen et al. 


2010; 


Fene: et al. 


2007) 



work with the joint observational data including the H[z) data are more stringent than 
the previous results for constraining the interaction term with other combined observations 
arising from the 182 Gold SNe la samples, the shift parameter of CMB given by the 3-year 
WMAP observat ions, the BAO measurement from SDSS and age estimates of 35 galaxies 



flFeng et al.ll2008[ ) 



4. Conclusions 

In this paper, we have examined, with the newly revised H{z) versus redshift data, as 
well as the baryonic acoustic oscillation (BAO) peak detected by large-scale correlation func- 
tion of luminous red galaxies from Sloan Digital Sky Survey (SDSS), the cosmic microwave 
background (CMB) observation and the 557 newly released Union2 SNe la data, to constrain 
two phenomenological interaction models for dark energy and dark matter, which are pro- 
posed as candidates to ease the coincidence problem of the concordance ACDM model. We 
find that, for the 7^ IDE and 7^ IDE models where 7^ and 7^ quantify the extent of interac- 
tion, although the H{z) data can not tightly constrain the model parameters, stringent con- 
straints can be obtained in combination with the CMB observation from the WMAP7 results, 
the BAO observation from the SDSS Data Release and the Union2 SNela set. In order to ex- 
amine the role of the H{z) data played in cosmological constraints, we also compare the SNe 
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Fig. 2. — The 68.3% and 95.4% confidence level contours for wx versus 7^ with the H{z) 
data (a); iy(2)+BA0+CMB (b); SNela+BAO+CMB (c); if(^)+SNeIa+BAO+CMB (d). 
The dashed lines represent 7^ = and wx = — 1- 
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la data in the same way and find the constraints with H{z)+BAO+CMB, SNe+BAO+CMB, 
and if(2;)+SNe+BA0+CMB combinations are consistent with each other. 



For the 7^ IDE model, we obtain wx = — I-IO^q;}^, 7^ 



which shows that the special case (7^ = 0, wx 



■ -0.013+^:^i^ with/f(2)+BA0+CMB, 
corresponding to the ACDM with no in- 



teraction) is excluded at la; whereas for SNe la+BAO+CMB and H{z)+SNe la+BAO+CMB, 



the best-fit values are 7„ 



-0.009: 



-0.013 
-0.012 



and 7„ 



-0.011^o;on, which indicates that 



the ACDM model is still included within la error region. For the 7^ IDE model, with 
H{z)+BAO+CMB, the best-fit values are (wx,7d) = (-1.24;°;}^, -0.040^°;°^^) ^nd the spe- 
cial case (7rf = 0,wx = —1) is excluded at la. However, with SNe Ia-|-BAO-|-CMB and 
H{z)+SNe la+BAO+CMB, the best-fit values are = -0.023l|]:||^^ and 7^ = -OMOtomj, 
respectively. The ACDM model is still included within la error regio n. These results 



Fene: et al. 


2008; 


Chen et al. 


2010) 



In summary, our results show that the concordance ACDM model still remains to be a 
good fit to the recent observational data. However, the interaction that the energy trans- 
ferring from dark matter to dark energy is slightly favoured over the interaction deliver- 
ing energy from d ark energy t o dar k matter, wh i ch is consistent with those obtained by 
Guo et al.l (|2007| ): iFeng et al.l (|2008[ ): IChen et al.l (|2010[ ). therefore, the coincidence prob- 
lem is quite severe in the two phenomenological scenarios. We are looking forwards to see 
whether these results may be changed with future observational data of high redshift SNe la 
data from SNAP etc (lAbraham et al.ll2004j ). more precise CMB data from the ESA Planck 



satell ite (iBalbi 



2010; 



Gao et al. 



(lAUen 



2OO2I: 



et al.l 



2007), as well as other complementary data, such as Gamma Ray Bursts 



200' 



Cao fcZhu 



2008 



20111). 



data (lSchaeferll2007l:lBasilakos fc Perivolaropoulosll2008t iLiang et al.ll2008l : iLiang. Xu fc Zhu 

I g a I OlOT) . the data of X-ray gas mass f raction in clusters 



Ettori et al.ll2009[ l and gravitational lensing data (1Zhulll998l ; ISereno 
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ii"(^)+BAO+CMB 


SNe la+BAO+CMB 


if(;2)+SNe la+BAO+CMB 


^m,0 
Wx 
Im 


The 7„ IDE Model 


^•^""^-0.030(-0.052) 
^+0.16(+0.23) 
-'-•-'-"-'-0.17(-0.24) 
n ni q+0.013(+0.019) 


^•^'^-0.029(-0.050) 
^r,+0.12(+0.20) 
-'^•^^-0.13(-0.20) 

U-UU»_Q Qj^2(-0.019) 


^„^.+0.030(+0.050) 
°-0.030(-0.050) 
1 np-+0.11(+0.17) 
i.UO_Q -^2(_o.l7) 

"^•"^-'■-'--0.011(-0.017) 


^m,0 
Wx 
Id 


The 7rf IDE Model 


-| r,^+0.17+0.28) 
-'-•^^-0.19(-0.28) 

u.U4:U_Q Q42(-0.062) 


"-'•^'^-0.030(-0.051) 
-, p|r)+0.09(+0.20) 
-'-•^^-0.09(-0.21) 


U-^Dy_0.030(-0.045) 
, -,p,+0.13(+0.19) 
-'-■-'-^-0.13(-0.19) 

^■^'^^-0.037(-0.056) 



Table 1: The best-fit value of parameters {^mo-, wx, 7m (id)} for the 7^ IDE, 7^^ IDE 
with 1-a and 2-cr uncertainty for i7(;2)+BAO+CMB, SNe la+BAO+CMB and H{z)+SNe 
la+BAO+CMB, respectively. 



